Properties of the strongly paired fermionic condensates 
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We study a gas of fermions undergoing a wide resonance s-wave BCS-BEC crossover, in the BEC 
regime at zero temperature. We calculate the chemical potential and the speed of sound of this 
Bose- Einstein-condensed gas, as well as the condensate depletion, in the low density approximation. 
We discuss how higher order terms in the low density expansion can be constructed. We demonstrate 
that the standard BCS-BEC gap equation is invalid in the BEC regime and is inconsistent with the 
results obtained here. The low density approximation we employ breaks down in the intermediate 
BCS-BEC crossover region. Hence our theory is unable to predict how the chemical potential and 
the speed of sound evolve once the interactions are tuned towards the BCS regime. As a part of 
our theory, we derive the well known result for the bosonic scattering length diagrammatically and 
check that there are no bound states of two bosons. 

PACS numbers: 03.75.Hh, 03.75.Ss 



I. INTRODUCTION 

In their seminal papers A. Leggett [3] and P. Nozieres 
and R. Schmitt-Rink 2] studied a system of spin- 1/2 
fermions with attractive short-ranged interactions in the 
singlet channel at low temperature. If the interac- 
tions are weak, the fermions form a Bardeen-Cooper- 
Schrieffer (BCS) superconductor. If the interactions are 
made progressively stronger, at some critical interaction 
strength a bound state of two fermions becomes possi- 
ble. These bound states are then bosons, which undergo 
Bose-Einstein condensation (BEC). Studies by these and 
many other authors demonstrated that as the interac- 
tion strength is increasing, the BCS superfluid smoothly 
evolves into the BEC superfluid, without undergoing any 
phase transition in between. 

Quite remarkably, Leggett, and Nozieres and Schmitt- 
Rink afterwards, found that the BCS gap equation, the 
mean-field equation which describes the properties of the 
BCS superconductor, evolves under the strengthening of 
the potential into the Schrodinger equation of a bound 
pair of fermions. This happens despite the fact that the 
gap equation is derived on the assumption of weak in- 
teractions between the fermions and should in principle 
break down as the interaction strength is increased. In- 
terestingly, this observation was first made as far back as 
in 1966 by V. N. Popov, but his contribution, at least 
within the context of the BCS-BEC crossover, was not 
noticed until later. Similarly, the BCS-BEC crossover 
was first studied in 1969 by D. M. Eagles [1| while the 
two channel model [5] was first studied by Yu. B. Rumer 

i- 

The fact that the gap equation describes the BCS su- 
perconductor well, combined with the fact that the gap 
equation is also well suited to describe the bound pairs 
of fermions on the opposite side of the crossover, and 
taking into account that there is no phase transition in 
between, allowed Nozieres and Schmitt-Rink to conclude 
that it may be used to interpolate between these two 
regimes, when the potential is neither too weak nor too 



strong and the fermions form an intermediate crossover 
BCS-BEC superfluid. 

Interest in this subject was recently revived when the 
BCS-BEC superfluid was obtained experimentally in the 
cold atomic systems where the interactions can be tuned 
with the help of Feshbach resonances d, 0, H, H, [13, EH ■ 
Thus for the first time theoretical predictions regarding 
the BCS-BEC crossover can be tested experimentally. 

Given the lack of justification for the gap equation in 
the BEC regime beyond predicting the binding energy of 
a pair of fermions, the question remained whether it can 
be used to predict other properties of the BEC conden- 
sate, such as the condensate depletion, excitation spec- 
trum and so on. 

In this paper, we show that the gap equation is actu- 
ally not valid on the BEC side of the crossover. Although 
correctly predicting the binding energy of the fermionic 
bound pairs, it fails to describe the interactions of these 
dimers properly. As a result, if taken literally, it incor- 
rectly predicts the physical properties of the BEC super- 
fluid. 

Indeed, it is generally believed in the literature P, Q 
that the transition between the BCS and the BEC phases 
occurs without a phase transition making these phases 
qualitatively the same phase. Thus on a qualitative level, 
it is possible to use the gap equation to obtain order of 
magnitude estimates of the parameters of the BEC phase. 
However, if one is interested in quantitative calculations, 
one needs to go beyond the gap equation. Barring the 
exact solution of the problem, which most likely is not 
possible here, the next best thing is to identify a small 
parameter and do the calculations as an expansion in 
powers of that parameter. The gap equation, when writ- 
ten on the BEC side of the crossover, implicitly assumes 
that the small parameter is the interaction strength be- 
tween the bosons. This assumption is invalid; it has 
been known for some time now that the Born approxima- 
tion fails to describe the interactions between the bosons 
properly p^ . 

In this paper, we identify as a small parameter the 



2 



so-called gas parameter of the Bose gas, and do all cal- 
culations as an expansion in powers of that parameter. 

The gap equation should be replaced by another equa- 
tion on the BEC side of the crossover. That equation, 
described below, cannot be derived exactly, but can only 
be obtained as an expansion in powers of the gas param- 
eter. In the lowest order approximation, it allows us to 
compute the chemical potential fit,, the speed of sound u 
in the superfluid, and the condensate depletion to be, at 
zero temperature. 
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Here 



Ub « 0.60 a 



(1) 



(2) 



is the scattering length of the bosonic dimers, whose ap- 
proximate relationship to the scattering length of a pair 
of fermions a has first been derived by Petrov, Salomon, 
and Shlyapnikov The coefficient 0.60 is approxi- 

mate in the sense of being determined only numerically, 
although the procedure to find it is in principle exact. 
Ub is the density of bosons, which in the lowest order 
approximation in density can be replaced by n/2, where 
n is the density of the fermions. riQ^b is the condensate 
density, nib — 2m is the mass of the bosons, where m 
is the fermion mass, fib is the chemical potential of the 
bosons, which can be related to the chemical potential fi 
of the fermions by demanding that /Jb — 2/i coincides with 
the binding energy released when the boson is formed. 
Throughout this paper we adopt units in which h = I. 

Our results Eq. ^ coincide with the behavior of the 
dilute interacting Bose gas [13]. We would like to em- 
phasize that they are not exact, but are obtained as the 
lowest order approximation in powers of the gas param- 
eter an 3. It is possible to use the techniques described 
in this paper and calculate the next order corrections to 
Eq. II]). They will no longer necessarily coincide with 
the next order corrections in the standard dilute Bose 
gas. We will not calculate them here, however, and limit 
our discussion with merely indicating how these higher 
order corrections can be obtained. 

As the BEC condensate is tuned towards the crossover 
BCS-BEC regime, a increases reaching infinity at the uni- 
tary point which lies somewhere in the intermediate BCS- 
BEC region. Thus the approximation used in this paper 
breaks down as the intermediate regime is approached. 
Our technique is unable to tell us anything about the 
crossover regime, and we will not attempt to study it in 
this paper. 

All calculations in this paper are done at zero temper- 
ature. We think that our techniques can also be adopted 
to finite temperatures, including studies of the critical 
superfluid transition temperature, either in the large N 



approximation or using some other technique. See also 
Ref. [1^] for more references on this subject. We will not 
attempt to do this in this paper. 

It is important to note that the interaction potential 
between fermions is characterized not only by the scat- 
tering length a, but also by the effective range vq. Un- 
like a which changes as the interaction potential is tuned 
through the crossover, rg remains roughly the same. If 
|ro| 3> ri~3 , then the fermions are said to be in the nar- 
row resonance regime. In that regime, the gap equation 
is a good approximation to the actual physics of the BCS- 
BEC crossover everywhere However, the wide res- 
onance crossover where |ro| ^ n~3, |ro| <C \a\ is more 
relevant for the current experiments, and it is the wide 
resonance crossover which received most attention in the 
literature. Everything we described above this paragraph 
pertains to the wide resonance regime. In order to be 
certain that we work in the wide resonance regime, we 
employ the one channel model which guarantees that tq 
is vanishingly small. 

An important part of the theory presented in this pa- 
per is the diagrammatic derivation of Eq. It has come 
to our attention that while our work was in progress, 
this was derived diagrammatically, independently of 
us, by I. V. Brodsky, M. Y. Kagan, A. V. Klaptsov, 
R. Combescot, and X. Leyronas, in Ref. [l^. The tech- 
niques these authors employed to arrive at Eq. ([2|) coin- 
cide with the ones discussed here. 

Finally, we note that a problem similar to the one stud- 
ied in this paper was examined by L. V. Keldysh and A. 
N. Kozlov in the context of excitons in Ref. [l^- Those 
authors, however, concentrated on the case of Coulomb 
interactions of interest for the physics of excitons, while 
we are interested in short ranged interactions, which have 
quite a different physics. In particular, the universal re- 
sults Eq. ([T]) do not hold for their system. Additionally, 
they, as well as Ref. [1], used an operator version of per- 
turbation theory which is hard to generalize beyond the 
lowest order. In contrast, our theory is based on a func- 
tional integral and is easily generalizable to arbitrary or- 
der. It is still remarkable that these authors understood 
the breakdown of the gap equation in the BEC regime 
long before this became an important issue in the BCS- 
BEC crossover literature. 

The rest of the paper is organized as follows. In section 
|TT]we introduce the formalism and discuss how the calcu- 
lations of interest to us can be performed. In section Hill 
we set up a perturbative expansion in powers of the gas 
parameter. In the next section HVl we perform the actual 
calculations and derive Eq. ^ . Section |V] discusses an 
alternative derivation of the results of section [TVl and the 
origin of the breakdown of the BCS-BEC gap equation. 
Finally, in section IVTl which is followed by conclusions, 
we compute the scattering amplitude between bosons and 
the bosonic scattering length Eq. ([2]) , a necessary step in 
the calculations performed in this paper. Some of the 
technical details of the calculations are discussed in the 
appendices. 
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II. FORMULATION OF THE PROBLEM 

Consider a gas of spin-1/2 fermions, interacting with 
some short ranged interaction in the s-wave channel. The 
most straightforward way to study this gas is by means 
of a functional integral 



Z : 



(3) 



where S*/ = 5*0 + S\n_t is the action consisting of the free 
part 

and the interaction part 

Sint = Ay d^x (it i/iii/ixViV'T- (5) 

This term describes interactions which happen at one 
point in space. We need to remember that by itself such 
an interaction is unphysical and has to be supplemented 
by a momentum cutoff A. A > to reflect that we choose 
an attractive interaction potential. 

As is standard in the treatment of superconductivity, 
we introduce the Hubbard-Stratonovich field A, which 
allows us to rewrite Eq. ([3]) in the equivalent way 



(6) 



where the action S is now given by S" = So + Shs j ^^nd 



•S'hs 



(Tx dt 



(7) 



The total action S is now quadratic in fermions, so the 
fermions can be integrated out. This results in the fol- 
lowing partition function 



Z = VAVA e^^^ , 



(8) 



where the effective action S'a is given by 



S'a = — « tr log 



-A 



2m 



/ dtd^x AA, 



(9) 



where w"*" ~ lo -\- iOsign oj. The action S'a and the func- 
tional integral Eq. ([8|) represent the starting point in our 
theory, although sometimes it is also advantageous to 
keep the fermions explicitly, as in Eq. © and Eq. ([7]). 

In the modern literature on Feshbach resonances, the 
system described in Eqs.dH]), ([4]), and ^ is often referred 
to as the one channel model. The interactions between 
fermions in this model are such that their effective range 



ro is vanishingly small (in the limit of a large momentum 
cutoff A). And indeed, the results described in this paper 
are applicable to the limit |ro| <C |ro| <C \a\ only, 

thus the one channel model is a natural starting point for 
our study. Alternatively, one could also study the two 
channel model The two channel model involves an 
additional parameter which allows control of the value of 
To independently of a. However, if we chose to work with 
the two channel model, to ensure the smallness of ro we 
would need to study the limit of infinitely strong coupling 
within that model. In that limit, the distinction between 
the two channel and one channel models disappear. This 
is explained in more detail in Appendix [X] Because of 
this equivalence, we will not study the two channel model 
any further in this paper. 

The standard approach at this stage is to find the ex- 
tremum of the action Sa with respect to the field A. 
Assuming that the extremum of the action is when the 
field A takes some constant value Aq, we find the BCS- 
BEC gap equation 



1 
A 



d^p 



1 



(27r)3 



\2m 



(10) 



The integral in the gap equation is up to momenta of the 
order of the cutoff A. It is advantageous to trade the 
interaction strength A for the fermionic scattering length 
a using the relation 



47r 
mA 
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TT 



(11) 



This relation is fairly standard. For completeness, we 
include its derivation in Appendix |B] Then the BCS- 
BEC gap equation acquires the form 



Ana 



d^p 
(2^ 



(p2_ 

y2m 



2m 

„2 



(12) 

The integral over momentum p is now convergent and 
can be extended to infinity. 

If a is negative and |a|n3 <c 1 (this corresponds to the 
weak attraction) , then we say that the gas of fermions is 
in the BCS state. In this regime, the gap equation has 
been investigated in detail in the literature devoted to 
superconductors. It was shown that the fluctuations of 
A about its mean fleld value Aq, which solves Eq. (|12p . 
are indeed small, and the physics of the superconductors 
can be captured by solving Eq. (fT^ . 

As the interaction strength is increased, a becomes 
more and more negative, reaching negative infinity at the 
interaction strength which corresponds to the threshold 
of bound state creation. For even stronger potentials, a 
becomes positive, first large and then it gets smaller. In 
this regime there is no justification for using Eq. ^ 
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Yet Ref. argued that in the "deep BEC" regime 
where an^ <C 1, Eq. ((T^) can be reinterpreted as the 
Schrodinger equation of a pair of fermions, with chemical 
potential playing the role of half of the energy. So the gap 
equation then correctly predicts the formation of bosonic 
dimers. Indeed, solving the gap equation in this regime 
by anticipating that /i <^ — Aq, we neglect Aq and find 
within this approximation 



(13) 



1/raa? is the binding energy of two fermions, as expected. 

Nevertheless there is no reason to expect that the gap 
equation holds beyond Eq. p^ . For example, one could 
expect that attempting to use it to compute the deviation 
of /i from Eq. due to nonzero Aq in Eq. would 
lead to an incorrect relation between the boson chemical 
potential /ib = 2/i+ l/(ma^) and Aq. On the other hand 
knowing the proper relation between /ifc and Aq would 
allow us to compute the physical properties of the BEC 
condensate, such as the speed of sound. So there is a need 
to develop a reliable technique which allows computation 
of these quantities. 

In this paper we investigate the properties of the sys- 
tem in the low density "deep BEC" regime, that is, in the 
regime where a is positive and small, so that arts ^ 1. 
To this end, we will calculate the normal and anomalous 
propagators of the bosonic A field using a diagrammatic 
expansion. In doing so, we will keep only the diagrams 
lowest in the density of the system. However, we will 
not assume that the interactions are weak, or in other 
words, that the interactions can be treated in the Born 
approximation. 

The fact that the Born approximation is not applicable 
to the bosonic dimers has been known for quite a while. 
Indeed, the dimer's scattering length in the Born approxi- 
mation is flb = 2a [Tsj , while the correct scattering length 
is given by Eq. Although first derived in Ref. |12 | 
by solving the four-body Schrodinger equation, it can 
also be derived by summing all the diagrams contribut- 
ing to the scattering exactly. In this paper we perform 
this summation and derive Eq. ^ diagrammatically in 
section IVH This allows us to recognize ab when the dia- 
grams contributing to it appear within the diagrams for 
the normal and anomalous bosonic propagators in our 
theory. 

Once the bosonic propagators are calculated, we will 
impose the condition that they have a pole at zero fre- 
quency and momentum, representing the sound mode of 
the condensate. This will give us a relation between the 
chemical potential /x, the scattering length a, and the ex- 
pectation value Aq of the field A, which will replace the 
gap equation Eq. (H^). In the dilute Bose gas literature, 
this is usually referred to as the Hugenholtz-Pines rela- 
tion . This equation replaces the gap equation in our 
theory. 

If instead of this proper procedure we had used the gap 
equation, we would have arrived at Eq. ([T]) with = 2a, 



which is the Born approximation to the dimer scattering 
length. Thus the gap equation fails because it cannot 
treat the interaction between dimers beyond Born ap- 
proximation. 

In order to be able to determine the unknown quanti- 
ties /I and Aq, one also needs another equation. Its role 
is usually played by the particle number equation, which 
states that the total particle density is equal to n. 



1 

V 



(14) 



where V is the volume of the system. In its most naive 
incarnation, we calculate the particle density by assum- 
ing that the field A does not fluctuate and is equal to its 
mean value Aq. Then we find 



(2^ 



1 - 



p 

2m 



2in 



^0 



(15) 



We are going to show that this equation is indeed correct 
in the lowest order approximation in density. However, if 
one wants to go beyond the lowest order approximation, 
Eq. (|15p acquires nontrivial corrections. 



III. DIAGRAMMATIC EXPANSION 

To set up a diagrammatic expansion, we first write 
down free propagators which follow from Eq. ([7]) or 
Eq. ([H]), taking into account that we work in the BEC 
regime where /i < 0. The propagator of the fermions is 
given by 



1 



2m 



(16) 



Note that since /i < 0, the propagator is retarded. 

The propagator of bosons is, naively, equal to —A. 
However, we also need to include in it the self energy 
correction due to fermionic loops. In other words, we 
need to expand the effective action up to quadratic 
order in A and the corresponding term in the expansion 
corrects the propagator. This expansion is performed in 
Appendix [B] We find, again taking into account that 
M < 0, 



1 



(17) 



21. 

4m 



2^ - iQ 



The bosonic propagator is also retarded, as any propa- 
gator of free bosons must be at zero temperature [l3|. 

The fact that all the free propagators are retarded cru- 
cially simplifies the diagrammatic expansion. Note that 
in the BCS regime where /Lt > this simplification would 
not have taken place. Similar effect happens in the more 
standard dilute Bose gas [l^ . 
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FIG. 1: The lowest order contribution in the low density ap- 
proximation to the anomalous fermion self energy, Ea. The 
dashed line is a condensate line. 

The diagrammatic technique must also involve the 
lines beginning or ending in the condensate. These are 
assigned the value Aq. Although Aq does not coincide 
with the density of particles in the condensate, it is pro- 
portional to it. In turn, the condensate density is always 
lower than the total density n. Thus the expansion in 
powers of the density n can be replaced by the expansion 
in powers of Aq. Since each line beginning or ending in 
the condensate is assigned the value Aq, the expansion 
in the density can also be understood as the expansion in 
numbers of lines beginning or ending in the condensate. 

Now we would like to construct self energy corrections 
to the propagators Go and I?o, to find the low density 
approximation to the exact fermionic normal and anoma- 
lous propagators G„ and Ga and exact bosonic normal 
and anomalous propagators Z3„ and Da- We need to 
know Gum to calculate the particle density Eq. pij) . We 
need to know Dn,a to relate the chemical potential ^ to 
other parameters in the theory. 

We denote E„_a the normal and anomalous self en- 
ergy corrections to Gn,a- Keeping with tradition p^ . 
we reserve the notation En and E20 for the normal and 
anomalous self energy terms in bosonic propagators -D„,q. 
We write E = E^^) + E^^) + . . where E^J) denotes the 
contribution to E from the diagrams with j lines begin- 
ning or ending in the condensate. 

The only diagram with one condensate line which con- 
tributes to Ea is shown on Fig[TJ We find 

E(i) = Ao. (18) 

The reason for the absence of any other diagrams at this 
order is the absence of vertex corrections in our theory 
at /i < when all the propagators are retarded. Thus it 
is a feature of the BEC regime only. 

There are no normal self energy terms with one exter- 
nal line, thus E^^' = 0. Computing the normal Green's 
function G„ with this one self energy correction is equiv- 
alent to computing the correlation functions in Eq. (|14p 
by assuming that A does not fluctuate and is equal to 
Aq. This gives Eq. (fT5)) . 

No diagrams with two condensate lines contribute to 
the anomalous self energy, and thus E^ =0. The dia- 
grams contributing to the normal self energy at this order 
are shown in Fig. [21 The square block shown on this fig- 
ure is equal to the sum of all diagrams which contribute 
to boson-fermion scattering. There is an infinite number 



FIG. 2: The next to leading order contribution in the low 
density approximation to the normal fermion self energy. 




FIG. 3: The normal (top) and anomalous (bottom) self energy 
contributions to the bosonic propagator. The wavy lines are 
bosonic propagators. 



of such diagrams, and they can all be summed by solv- 
ing the Lippmann-Schwinger integral equation. This will 
be discussed in more details in section I VII For now we 
denote the result of summation as t^^. Then the contri- 
bution to the normal self energy at this order is given 

by 

El2)=A^t^^ (19) 

If needed, diagrams with three or higher number of 
external lines can also be constructed. 

Turning to the bosonic propagator, we find that there 
are no diagrams contributing to En and E20 at the order 
of Aq. At the order of Aq, the self energy diagrams 
are shown on Fig [3l The square box denotes the sum 
of all the diagrams which contribute to the T-matrix of 
scattering of a boson off another boson. Just as those 
for t^^ , these diagrams can also be summed up using the 
appropriate Lippmann-Schwinger equation. We denote 
the result of summation t^^. t^^ is calculated in section 
Thus we find 

E(^,U2Agibb, 4o^=Ag^^^ (20) 

The coefficient 2 appears as the appropriate combinato- 
rial factor. 

It is also possible to consider diagrams higher order in 
Ao which contribute to the bosonic normal and anoma- 
lous self energy. One such diagram, which contributes to 
E^^\ is shown on Fig (4] There are infinitely many simi- 
lar diagrams contributing to E^^"* and E20''. We will not 
discuss these any further. 
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FIG. 4: A possible higher order contribution to the bosonic 
self energy. 



IV. CALCULATION OF THE SPEED OF 
SOUND AND THE CONDENSATE DEPLETION 

Once the self energy terms are known, the procedure 
for calculating the chemical potential, the speed of sound, 
and the condensate depletion is fairly standard. 

First of all, we calculate the normal and anomalous 
bosonic Green's functions Dn,a, given En and E20. The 
calculation involves solving the appropriate Dyson equa- 
tion ^3j] and the result is 



S20(P,^) 

D{p,iu) ' 



D{p,uj) 



(21) 



where 

D{p,iu) 



- [S20 ip, Lof + {D^\p, Lo)-Enip,to))x 
X {D^\-p, -c^)-En(-p,-c^)). 

(22) 



The Hugenholtz-Pines relation then takes the form 

15(0,0) = 0. (23) 

In the lowest order in density, the self energy terms En 
and E20 vanish. Then the Hugenholtz-Pines relation 
reads 



^o"'(0,0) = 0. 



(24) 



Recalling the definition of Do, Eq. (fT7|) . we immediately 
find fi = — l/(2ma^), the same relation as the one which 
follows from the gap equation, Eq. 

In the first nonvanishing order in density, the self en- 
ergy terms are given by Eq. ((20|l . We are interested in 
the self energy at zero momentum and frequency, thus 
t^^ also needs to be computed at zero momentum and 
frequency. In section fVll we show how to calculate t^^ in 
vacuum by solving the appropriate Lippmann-Schwinger 
equation. We, however, need t^^ at a finite chemical po- 
tential /X. Fortunately, if 



l/{2ma^)\ < \n\ 



(25) 



then t^^ at a finite /i is approximately given by its ex- 
pression computed in vacuum at zero momentum and 



frequency. In other words, it is proportional to af,, where 
Of, is the bosonic scattering length in vacuum, Eq. 
The Lippmann-Schwinger equation can only be solved 
numerically, so we only know the numerical value of ai,. 

More precisely, t^^ at zero momentum and frequency, 
and at chemical potential satisfying Eq. ((25)) is related to 
the scattering length ai, as follows 



bb 



27rab 



(26) 



where Z is the residue of the bosonic propagator Dq at 
its pole, or 



(27) 



We expect Eq. (^5]) to hold as we expect fi to deviate 
only slightly from its zero density value Eq. (|13[) in the 
low density approximation employed here. 

Armed with these relations, we introduce a bosonic 
chemical potential 

1 



= 2^ 

flLU- 

and solve the Hugenholtz-Pines relation in the next order 
in density 



fib = 



(28) 



There are actually two solutions of the Hugenholtz-Pines 
relation, but one of them is well known to be unphysical 

El. 

We remark that if we had instead decided to solve the 
gap equation Eq. ([12]) by expanding its square root in 
powers of Ag, we would have obtained fib — AQa^m/2, 
as if the bosonic scattering length were ai, = 2a. This is 
manifestly incorrect. Thus we see that the gap equation 
fails in the EEC regime. In fact, it is easy to see exactly 
where the gap equation breaks down, and what comes in 
its place. This analysis is performed in section [V] 

Now we use the particle number equation Eq. (jl5p to 
relate Aq to the particle density n. Within the low- 
est approximation in density, we can substitute /i — 
-l/(2ma2) to find 



am^ Aq 
47r 



(29) 



Using the next order approximation for the chemical po- 
tential would only contribute to the further expansion of 

Eq. in powers of Aq, where, however, we would also 

(2) 

need to take into account the terms Ea,ri. Fortunately 
this is not needed for the calculations in the lowest order 
presented here. 

Combining Eq. ^ and Eq. ^ we find 



m rub ' 



(30) 
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where rub = 2m is the mass of the bosonic dimers and 
nt = n/2 is their density. This is the first of the resuhs 
advertised in the beginning of the paper in Eq. ([T|). 
To find the speed of sound, we impose the condition 



D{p,u) = 0. 



(31) 



This gives us a relation between lu and p. For small lu 
and p, with the help of Eq. ([29| it reduces to 



(32) 



2 2 47rnbafe 
w =p —2—. 



Comparing with = u^p^, where u is the speed of 
sound, gives the second result in Eq. ([1]). 

Finally, we calculate the condensate depletion. The 
density of bosons not in the condensate can be calculated 
as 



Srib = — hm 



dio (Pp 



Dnip,Lu) e 



(33) 



The calculation of the integral in Eq. ([33|) demands know- 
ing the full frequency and momentum dependence of the 
self energy terms Sn and E2o. This in principle can 
only be found numerically, since the knowledge of full 
frequency and momentum dependence of t^^ is required. 
However, in the limit of small density where self energy 
terms are small, the relevant frequency and momentum 
which contribute to the integral in Eq. (|33p are in the 
range uj ~ Aq and p^ /{Am) ~ Aq. Thus we expand 
^ in powers of A§ ^ \uj +p^/(2m;,) — <C l/{ma?) 
to find the approximate expression for the free bosonic 
propagator Eq. ([TT]) 



Z 



1 



to ■ 



4m 



(34) 



At the same time ZS20 coincides with its expression for 
the dilute interacting Bose gas. With this substitution, 
Dn{p, oj) is equal to that of the dilute Bose gas, with the 
scattering length at. Thus the standard formula for the 
condensate depletion [l^ holds, and it reads 



rib al 



(35) 



which is the final result in Eq. ([T]). We emphasize that 
this is also an expansion in powers of na^, and that next 
order terms will likely depend on the functional depen- 
dence of the self energy terms on the frequency and mo- 
mentum. 

Let us summarize what one needs to do to go beyond 
the lowest order approximation discussed above. The di- 
agrams for the higher order corrections to Sn and S20J 
such as the one shown on Fig. |4l have to be evaluated. 
There will be an infinite number of such diagrams, and 
they can possibly be summed up using a numerical tech- 
nique similar to the one employed here to calculate t^^. 



One also needs to calculate higher order corrections to 
Sq and S„. In particular, at the order of Aq, t^^ has to 
be evaluated. Because of the nonzero chemical potential 
present in the external lines of the diagrams shown on 
Fig. [2I t^^ will not reduce simply to the scattering length 

abf, and so extra work will be needed to find the value 
(2) 

of S„ . As pointed out by Hugenholtz and Pines in Ref. 

, at each order one also needs to correct the propaga- 
tors such that a gapless sound mode is produced at this 
order. 

Once the self energy terms are found in the desired 
order, one needs to solve the Hugenholtz-Pines relation 
in the next order of Aq. This also involves taking into 
account the dependence of£°^ on ^b, which was neglected 
here due to fib being zero up to terms of the order of Ag. 
This allows to find the chemical potential. Given the 
chemical potential, one can also find corrections to the 
spectrum by looking at the poles of the propagator. For 
that, the expansion of t^^ in powers of momentum will 
also be needed. This expansion is also required in order 
to find corrections to the condensate depletion, given at 
the lowest order in Eq. (pSj). 



V. THE ORIGIN OF THE FAILURE OF THE 
BCS-BEC GAP EQUATION 

An alternative method of deriving the results of the 
previous section is presented here. A convenient way to 
derive the BCS-BEC gap equation is by shifting A by 
Aq in Eq. ([6]) and differentiating the resulting Z with 
respect to Aq, assuming that Aq gives an extremum of 
logZ. This gives 



^ = (V'iV'T)' 



(36) 



To derive the naive gap equation Eq. ^12\\ one usually 
evaluates the right hand side of Eq. ([36| at constant A = 
Aq. While this procedure works in the BCS regime, it 
actually breaks down in the opposite BEC regime. 

Indeed, on the BEC side of the crossover, it is more 
appropriate to expand everything in powers of Aq. So we 
instead compute the right hand side of Eq. ((36)) in terms 
of diagrams with a fixed number of condensate lines. 

The simplest diagram contributing to Eq. (|36p is shown 
on Fig. [5l There are two condensate lines shown on that 
diagram. The line on the left is not an actual line. Rather 
it is drawn simply for convenience, and the vertex where 
it is attached to the diagram represents V'iV'T whose cor- 
relation function we are computing. The line on the right, 
however, is the actual condensate line. Thus this diagram 
has one condensate line and is proportional to Aq. 

If we use this diagram for the right hand side of 
Eq. we derive Eq. with the substitution Aq = 0. 
Thus this reproduces the lowest order gap equation, and 
gives Eq. (fTS]) for fi. 

This is the only possible diagram with one condensate 
line. There are no diagrams with two condensate lines 






FIG. 5: The simplest diagram contributing to the gap equa- 
tion Eq. pel). 



FIG. 8: Diagrams contributing to the f-matrix of fermion- 
boson scattering. 



either. Once we allow for three external condensate lines, 
we find one possible diagram shown on Fig. [S) 




FIG. 6: The simplest diagram contributing to the gap equa- 
tion at order Aq. 

This diagram, when taken into account, gives the usual 
gap equation Eq. (|12p where the square root has been 
expanded up to linear order with respect to Aq. We 
know that solving Eq. at this order gives a?, = 2a 
And indeed, the diagram shown on Fig.[6]is nothing 
but the simplest process contributing to the boson-boson 
T-matrix, which is the very first diagram on Fig. 1111 In 
other words, it is the Born approximation for the boson- 
boson scattering. 

Wc know that the Born approximation breaks down 
in the BEG regime. Instead, at this order we need to 
take all the diagrams with three external lines, which are 
equal to t^^, as shown on Fig. [71 This gives the correct 



■>- 




->- 



FIG. 7: The sum of all diagrams contributing to the gap 
equation at order Aq. 



bound states of fermions. All calculations proceed in 
vacuum so we set ii — for this entire section. As ar- 
gued in section IIVI these vacuum scattering processes 
are precisely those needed to calculate the normal and 
anomalous fermion self energies at next to leading order 
in the low density expansion and the lowest order boson 
normal and anomalous self energies, respectively. 

The scattering t-matrix consists of all diagrams with 
incoming and outgoing lines corresponding to the scatter- 
ing particles. Since we work in vacuum, there is no hole 
propagation, or all the propagators are retarded. These 
observations greatly simplify the possible diagrams. In 
fact, had any of the propagators been advanced, there 
would be no hope of summing the diagrams contributing 
to the scattering processes exactly. Fortunately, when all 
the propagators are retarded, we are able to sum all the 
diagrams using Lippmann-Schwinger type integral equa- 
tions. We emphasize that all the calculations proceed 
in the regime of vanishing effective range tq, so the re- 
sults of this sections are valid only for potentials where 
|a| » ko|. 

Wc first consider scattering of a fermion and a bound 
state of two fermions, a bosonic dimer. This process was 
first solved in the zero effective range approximation by 
Skorniakov and Ter-Martirosian p3| in the related prob- 
lem of neutron-deuteron scattering. We include here a 
somewhat detailed description of this scattering process 
in part because it is a useful exercise before consider- 
ing the more complicated case of scattering between two 
bosons. 



p, p^^]ilAm-h 



k. k'lAm-h 



— ? — 










— > — 



£ k^lAm-b 




—p,—p^+k /2m 



k-l2m pi p^+lcl4m-b 



gap equation, with 2a being replaced by 0.60 a. 

To summarize, instead of using the naive gap equation 
Eq. ([T2I) which does not take into account fiuctuations 
around Aq, in the BEG regime we use Eq. (|36p . expand- 
ing the right hand side in powers of Aq. This procedure 
goes beyond the Born approximation and reproduces the 
correct results in the BEG regime. 




p, p^^+k lAm-b 



FIG. 9: The integral equation for the t-matrix of fermion- 
boson scattering. 



VI. SCATTERING OF PARTICLES AND 
BOUND STATES 

In this section we consider the scattering of a fermion 
by a bound state of fermions and scattering between 



The first few diagrams contributing to the t-matrix are 
shown in Fig. |8l where the external lines are not part of 
the t-matrix. Even though all the propagators are re- 
tarded, there is still an infinite number of diagrams con- 
taining retarded propagators only, and they all contribute 
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at the same order. Indeed, from these first examples it 
is clear that any such contribution to the i-matrix will 
have n > boson lines of order a/m, see Eq. (fT7| . 2n+l 
fermion lines of order ~ ma^ and n integrations of 



order Ep'^ 



Thus each diagram contributing 



to the f-matrix is of order ma . 

There is no hope of computing these diagrams one by 
one and summing them up. Instead, we derive an integral 
equation which their sum satisfies. Let the incoming bo- 
son and fermion have on-shell 4- momenta (fc, k"^ /Am — b) 
and (— fc,fc^/2m), respectively, and let the outgoing bo- 
son and fermion have 4-momenta {p,pQ + k'^/ Am — h) and 
{—p,—po + fc^/2m), respectively. Here b is the binding 

—b. The t-matrix with 



energy of the dimer, —1/ (ma" 
these kinematics is denotec 
satisfies the integral equation 



these kinematics is denoted t^^{p,po). Then the t- matrix 



q(.P,Po) = -Goik+p,-kyAm-b + pa) 



xGo(-g, qo)Go{p + q,Po + qo - b), (37) 

Zm Am 

which is depicted in Fig. [9] The minus sign in front of 
the right hand side of the equation is due to fcrmionic 
anticommutations. The reason for letting outgoing mo- 
menta be off-shell is that this makes it possible to solve 
the integral equation. The on-shell point has |fc| = \p\ and 
Po = 0. It is possible to integrate out the loop energy on 
the right hand side in Eq. ([57|) by noting that ^^((f, go) 
is analytic in the upper half plane in go which can be 
easily seen by looking at the diagrams contributing to 
the t-matrix. This integration sets qa = k"^ /2m — q^ /2m 
and we thus set po = k^ /2m — p^ /2m, to have the same 
dependence in the t-matrix on both sides in the equation. 
Define tf{p,po = P/2m- p^/2m) = tf{p). 

At low energies scattering is dominated by s-wave scat- 
tering, thus Eq. (|37p can be averaged first over directions 
of k and then over directions of p. The angular average 
of t^{p) is denoted t^{p)- Then the integral equation for 
the t-matrix of fermion-boson scattering becomes 



tl\p) 



m p 
■ m — 



- pk -\- k"^ — mE 



2pk 
1 



pk + k'^ — mE 



- In ^ 



qp + p^ — mE 



qp + p^ — mE 



- ^^J-E + iq'^/Ar 



(38) 



Here, E = 3fc^/4m — & is the total energy which is as- 
sumed to be less than zero. To calculate the scatter- 
ing amplitude each external bosonic leg has to be renor- 
malized by the square root of the residue of the pole of 
the bosonic propagator at the energy of the bound state, 
a) , compare with Eq. (|27p . The external 
fermionic legs are free propagators and do not have to be 



renormalized. The scattering amplitude is evaluated on 
shell and is 



T'^'{k)=Ztf{k) 



hi. 



(39) 



The relationship between the scattering length and the 
scattering amplitude is 



37r 



T-(0) = -a,/. 
m 

Solving Eq. pS)) it is found that abf ^ 1.18a. 



k, k'lAm-b p, p^+lcl^^m-b k, krlAm-b p,p^^+k/4m-b 



(40) 
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-k. k'lAm-b -p,-p^^+k^/4m-h -k, k^Am-b -p.-p^^+k lAm-b 

k,k'IAm-b q, q^^+k~/4m-b p, p^^+k~/4in-b 






r 




-q^^k'l4m-b -p, -p^^+k''l4m-b 



t/4m-b 

FIG. 10: The integral equation for the t-matrix of boson- 
boson scattering. F is the sum of all two boson irreducible 
diagrams. 

We now turn to the scattering between two bosons 
where the bosons are bound states of two distinguishable 
fermions. This process was first solved by Petrov et al [l^ 
by solving the quantum mechanical 4-body problem. As 
we already mentioned, recently, while the work reported 
here was in progress, the problem was solved in Ref. [lB| 
by using the same diagrammatic approach as we employ 
below. 

As in the case of fermion-boson scattering there are no 
condensate lines internally in the diagrams and no hole 
propagation. 

Consider the different diagrams which contribute to 
the t-matrix of boson-boson scattering. The t-matrix 
consists of all possible diagrams with two incoming and 
two outgoing bosons, the external legs not included. The 
crucial point is that all of these diagrams are of the same 
order, namely they are all proportional to m^a'^ which 
we will show below. This means that the diagrams con- 
tributing to the t-matrix do not form a perturbation se- 
ries and it is not adequate to keep only the Born approx- 
imation to the t-matrix. Instead an infinite number of 
diagrams must be taken into account. The summation 
of all diagrams contributing to the t-matrix may be per- 
formed by using an integral equation. In particular, the 
t-matrix may be built from a series of two boson irre- 
ducible diagrams, which results in an integral equation 
as shown in Fig. 1101 Two boson irreducible diagrams are 
understood to be diagrams which cannot be cut in half 
by cutting two boson lines only. 
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r 






FIG. 11: The simplest diagrams contributing to the two boson irreducible diagram, F. 
side corresponds to F*"' defined in Eq. (|Clip . 



The first diagram on the right hand 



Let the kinematics be as depicted on the figure with the 
incoming energies and momenta chosen on-shell and the 
outgoing energies and momenta ofi-shell to allow for the 
solution of the integral equation. The on-shcU condition 
is \p\ = \k\ and po = 0- The i-matrix with kinematics 
as shown on the left hand side of Fig. [10] is denoted 
by £^{'p,pq). The two boson irreducible diagram with 
(±g, E/2± qo) incoming 4- momenta and (±p, E/2ztpo) 
outgoing 4- momenta is given by Ffc((f, qo]p,po). E is the 
total energy, E ~ — 26 + The integral equation 

corresponding to these kinematics is 



tf{p,Po) 



+i 



Ffe(A:,0;p,po 



Doiq, E/2 + qo)Doi-q, E/2 - qo) 



(41) 



Thus far the integral equation does not depend on 
whether the constituents of the dimers are bosons or 
fermions. Since the scattering is dominated by s-wave 
scattering we average Eq. (PT|) over directions of k, then 
directions of p and finally the integration over directions 
of q can be trivially performed. Let the corresponding 
angular averages be t\^{p,Po) and Tkiq, qo;p,Po) and let 
k = \k\, p = \p\, and g = Then the integral equation 
becomes 



t't{p,po) = rk(k,0;p,po) + - 



^k{q,qo;P,Pa)tkiq,qo)q dq dqo 



1 + V(9V4 - E/2r - q^ - V2jqy4 - E/2 + ^{q^A- E/2)^ - q^ 



(42) 



r 



which contains no remaining angular dependence and 
where the product of bosonic propagators has been writ- 
ten explicitly using Eq. ([TT]). Since all diagrams are of 
the same order of magnitude, it is convenient to let all 
momenta be measured in units of and energies in 
units of 6 = l/(ma^). t^^ and F arc then measured in 
units of m^a^. For simplicity this is the case in Eq. (f42|) 
and will be the case in the remainder of this section. 

The simplest two boson irreducible diagrams are shown 
in Fig. [11] where the external lines are not a part of F. 
It is clear that there is an infinite number of such dia- 
grams. Letting n count the number of bosonic propaga- 
tors in the diagram contributing to F, any such diagram 
has A + 2n fermionic propagators, each of order ma^, n 
bosonic propagators of order a/m and n + I integrations 
of order m~^a~^. Thus each of these diagrams is of order 
m^a'^. To see that t^^ is also of order m^a'^ note that any 
diagram contributing to t^^ will contain n + 1 factors of 
F, 2n bosonic propagators, and n integrations, resulting 
in any diagram contributing to t^^ being proportional to 



Letting the first diagram on the left in Fig. [TJequal the 
t-matrix corresponds to the Born approximation which 
gives Of, = 2a. Using this same diagram for F reproduces 
the result of Ref. 21], a?, = 0.78a. However in both 
of these methods an infinite number of diagrams of the 
same order as those considered is ignored. 

The sum of two boson irreducible diagrams may again 
be obtained through solving an integral equation. The in- 
tegral equation to be solved is depicted in Fig. [T^ where 
S denotes all two boson irreducible diagrams with two 
incoming bosons, one outgoing boson and two outgoing 
fermions. The two boson irreducible diagram F may then 
be obtained by tying together the two external fermion 
lines. Furthermore, since the two boson irreducible dia- 
gram F does not depend on angles, we choose to let S be 
averaged over the direction of the incoming momentum. 
Define a four vector K = (0, £'/2) which describes the 
incoming energy carried by a boson. Then the integral 
equation corresponding to the kinematics shown on the 
figure is 
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FIG. 12; Integral equation for the two boson irreducible diagram with two bosons coming in and one boson and two fermions 
going out. 
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The minus signs in front of both the first order terms and 
the higher order terms are a result of anti-commuting 
fermions. The integral over dVL^ is the averaging over 
angles of incoming momentum. F is then found from S 
by connecting the outgoing fermion lines and integrating 
over the loop momentum. The precise relation between 
r and S is 

r.(g;p) = \ j 0^Go{K/2 + q')GoiK/2+p-q') 
xEkiq,q',p-q'). (44) 

It is not essential to include both lowest order contribu- 
tions to S in Eq. (|43p . the only difference if only one of 
these were included is that then the factor 1/2 in Eq. ([44]) 
should be removed. However, the symmetric structure of 
Eq. (1131) means that 

2fc(g,Pi,P2) = Sfc(g,p2,Pi) (45) 

which will be useful in the solution of these equations. 
How we solve the set of equations (li^ . pS)) . and is 
described in detail in Appendix [Cl 

As in the case of fermion-boson scattering, to calculate 
the scattering amplitude each external bosonic leg has to 
be renormalized by ^/Z = Q-k / {rm? a) . The scattering 
amplitude is then evaluated on shell 

T^^(k) = Z^t^t}'{k,G). (46) 

with k = \k\. Here units are restored to the t-matrix. 
The scattering amplitude is related to the boson-boson 
scattering length by 

Tbb(o) ^ ^a,. (47) 
m 

By solving Eq. (|^^ using the method described above 
and in Appendix ICl w e find « 0.60a in complete agree- 
ment with Refs. |12|. [T6j. 

Using the above formalism it may also be checked that 
there are no bound states of a pair of bosons. If such 



bound states exist, the interaction between the bosons in 
this theory can become effectively attractive, rendering 
the bosonic gas unstable. A bound state of two bosons 
would correspond to a pole in the scattering amplitude, 
that is, to a solution of the homogenous version of the 
integral equation ((4T|) for the t-matrix. Varying the total 
energy, that is /c^ for fc^ < 0, we do not find any solutions 
to the homogenous equation and thus do not find any 
bound states. 

It should be mentioned that the interaction between 
bosons does not have to be renormalized, in contrast to 
the interaction between bound states of bosons [l^l . In 
the latter case, the interaction between a boson and a 
dimer of bosons depends on the ultraviolet cutoff and a 
possible solution to this difficulty is to add a three body 
force counterterm to the theory and use an additional 
input, such as an experimentally measured three body 
scattering length [2^. However, this difficulty is absent 
in the case of scattering between a fermion and a dimer 
of fermions and also in the problem of scattering between 
two dimers of fermions. 



VII. CONCLUSION 

In this paper we discussed the low density expansion 
of the BCS-BEC condensate in the BEC regime. We 
found the dispersion of its Bogoliubov modes (the speed 
of sound), its chemical potential, and the condensate de- 
pletion in the lowest order approximation in powers of 
the gas parameter ana. Notice that the gas parameter 
increases as the system is tuned towards the BCS-BEC 
crossover regime, eventually reaching infinity at the so- 
called unitary point lying in between the BCS and BEC 
condensate. So the theory developed here works in the 
BEC regime only and breaks down in the crossover area. 

Throughout the paper we emphasized that the BCS- 
BEC gap equation actually breaks down in the BEC 
regime. We would like to remark further on the origin 
of this breakdown. The gap equation is derived by min- 
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imizing the effective action of tfie condensate, given by 
Eq. This is correct only if the fluctuations about 

this minimum are smafl. If they are not small, the fluc- 
tuations must be taken into account. In that case, min- 
imization of the effective action should be replaced by 
minimization of the effective potential, which takes into 
account fluctuations, as discussed in any textbook on 
field theory. In practice, it is possible to replace the 
evaluation of the minimum of the effective potential by 
the Hugenholtz-Pines relation, as is done here. The ad- 
vantage of this second procedure is in the fact that it 
automatically gives the excitation spectrum in addition 
to the chemical potential. The two techniques are how- 
ever equivalent. Most importantly, the calculations per- 
formed in this paper demonstrate that in including the 
fluctuations into the effective potential, one needs to go 
well beyond the Gaussian fluctuations approximation. In 
fact, terms up to infinite order (in terms of the number 
of loops) have to be summed up, all of them being of the 
same order in the gas parameter. Fortunately, this is pos- 
sible to do, and this is what has been accomplished in this 
paper. Alternatively, the gap equation corrected by fiuc- 
tuations was calculated in section |Vl We demonstrated 
how the naive gap equation breaks down due to fluctu- 
ations, and how the fluctuations effectively replace the 
Born approximation with the full boson-boson scattering 
amplitude. If we wanted, we could follow the procedure 
outlined in section|V]as an alternative to the Hugenholtz- 
Pines relation. 



The same argument also demonstrates that the Gross- 
Pitaevskii equation of the condensate should follow not 
from the effective action, but from the effective potential. 
In other words, the true bosonic scattering length Eq. ^ 
should be used in its quartic term, as opposed to its Born 
approximation value 2afc. The Gross-Pitaevskii equation 
will of course be valid only at length scales much bigger 
than a. 



All calculations in this paper are done in the lowest ap- 
proximation in density, which significantly simplifies the 
work needed to be done. Finding higher order corrections 
to Eq. ID) is one possible direction of further work along 
the lines discussed in this paper. Another possible di- 
rection is to study the condensate at finite temperature. 
This should probably be done in the large N approxima- 
tion such as the one used in Ref. [l^]- Generalizing the 
techniques of our paper to finite temperature should be 
a promising direction of further research. 



On the other hand, we do not expect that these tech- 
niques will help to shed light on the BCS-BEC crossover 
regime, especially at the unitary point. The small pa- 
rameter utilized here becomes infinity there. The unitary 
regime can only be understood numerically psj . barring 
an invention of an exact solution. 



Acknowledgments 

We thank J. Shepard, L. Radzihovsky, D. Sheehy, and 
A. Lamacraft for useful discussions. This work was sup- 
ported by the NSF grant DMR-0449521. J. L. also wishes 
to thank the Danish Research Agency for support. 



APPENDIX A: THE TWO CHANNEL MODEL 

The two channel model is defined by its functional in- 
tegral 



Z 



where the action Stc is given by ^tc = Sq + Sob + 
Here Sq is the free fermion action given in Eq. (|4]). 
is the free boson action given by 



(Al) 

- Sab- 
Sob 



s 



Ob 



d-^x dt b 



d 



1 92 



dt Am d:x? 



2^l-eA b, (A2) 



where eo is the detuning, the parameter controlled by 
the magnetic field in the Feshbach resonance setup. Sab 
is the interaction term 



Sab 



g / d^x dt {bip-^ipi +bipitp^) 



(A3) 



The scattering of fermions in the two channel model is 
characterized by the scattering length a and effective 
range ro given by [26l | 



[^0 - ^ ) 



ro 



An 



Sir 



(A4) 



In order to ensure that we work in the wide resonance 
regime |ro| <C n~3, |ro| <C |a|, in this paper we need to 
take the limit g 00, while simultaneously adjusting cq 
so that a remains finite. Introducing the notation A = 
gb, A — gb, we find that g disappears from Sab, while 
S'ob, in the large g limit, becomes 



Snb — ' 



m mA\ 



d-^xdtAA. 



(A5) 



We recognize the combination of parameters in brackets 
as — 1/A, thanks to Eq. (fTTj) . Thus the two channel model 
reduces to the one channel model as given by Eq. ((4]) and 
Eq. Q. 

One important lesson which follows from this discus- 
sion is that in the wide resonance regime, it is not justified 
to use g as a small parameter and construct perturbative 
expansion in its powers. Indeed, g is not only not small, 
it should actually be taken to infinity. 
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FIG. 13: The bubble diagrams of fermion scattering. 



APPENDIX B: FERMIONIC SCATTERING AND 
RENORMALIZED BOSONIC PROPAGATOR 

The scattering of fermions governed by the action Eq. 
(Ill) and interacting via the short range interaction Eq. ([5|) 
can be calculated by summing up the diagrams depicted 
in Fig. 1131 The calculation proceeds in vacuum, so the 
chemical potential /.t is set to zero everywhere. This sum 
gives the fermion T- matrix, , and the scattering length 
is proportional to the T-matrix at zero momentum, 



T^{0). 



(Bl) 



The bubble diagrams form a geometric series which can 
be summed to give 



T^(0) = 



-A 



1 + An(o) 



(B2) 



where IT is a bubble as shown in Fig. [13] evaluated at 
zero momentum. This bubble is given by 



n(0) = i 



Go{p)Go{-p) 



mA 
'2^' 



(B3) 



Here, A is the momentum cut-off. This shows Eq. dTTI) . 
namely that 



47r 2A 
mA TT 



(B4) 




The fermion loop can be calculated to give 



U{p) 



j^;^Gq{p + q)Go{-q) 
,3/2 



mA III ' I — 
27r2 47r 



392/4m - 2/^ - iO. (B6) 



Inserting this result in Eq. (jB5[) . using Eq. (|lip to write 
the momentum cut-off in terms of a, finally gives 



47r 



m _i 
a ^ — 



Am I 



(B7) 



Note that it is crucial that /x < for the integral in 
Eq. (jB6|) to be calculated the way it is. Otherwise, the 
fermionic propagator is no longer retarded due to hole 
propagation. Thus Eq. (jB7P applies only in vacuum or 
in the BEG regime where /x < 0. 



APPENDIX C: THE INTEGRAL EQUATION 
FOR BOSON-BOSON SCATTERING 



In this appendix we describe the solution of Eqs. 
B4)) . describing scattering of two bound states, each con- 
sisting of two fermions. We demonstrate how to inte- 
grate out the loop energies in constructing the sum of 
two boson irreducible diagrams, F, and we write the cor- 
responding equations in a more convenient way for nu- 
merical studies. We also describe the numerical methods 
used. For simplicity, in this appendix all momenta are 
measured in units of a^^, all energies in units of 1/ {ma?), 
and and F in units of 'm?a?' . 



in 



rv 



FIG. 14: The bosonic propagator renormalized by fermion 
loops. The thin wavy lines are unrenormalized boson propa- 
gators. 

The calculation of the bosonic propagator corrected 
by fermionic loops proceeds similarly. Having in mind 
its applications in this paper, we do this calculation at 
a finite chemical potential ^ < 0. Expanding the action 
5a in Eq. ([9]) results in a geometric series as shown in 
Fig. 1141 which can be summed to give 



FIG. 15: The original and rotated contours of frequency inte- 
gration used in Eq. (|42p . The crosses correspond to poles of 
the bosonic propagators which approach the imaginary axis 
as k,q ^ 0. The dots and black lines towards ±oo are branch 
cuts of the bosonic propagators and the grey boxes are pos- 
sible non-analytic structure of V{q,qo;p,po). These do not 
come closer to the imaginary axis than ±6. 

To construct Tk{q, qo;p,pQ) it is convenient first to ro- 
tate the external energies po and qo onto the imaginary 
axis which can be done without crossing any poles. Both 
applications we have in mind, calculation of the scat- 
tering length and search for bound states, require the 
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t-matrix t\^{p,po) evaluated on shell and thus we are 
only interested in the t-matrix evaluated at pq — which 
is still on the contour. That it is possible to rotate the 
external energies is not immediately obvious. Consider 
Eq. (|42p which relates the i-matrix to F. The two 
bosonic propagators have branch cuts in qo starting at 
go = ±(1 + — fc^/4) and going towards ±oo, re- 
spectively. l + g^/4 — /c^/4 > since to calculate the 
scattering length we let fc = and to search for bound 
states we let fc^ < 0. The bosonic propagators also have 
poles at qq = ±((7^/4 — fc^/4 — iO) where the infinitesi- 
mal is determined from the requirement that Eq. (jl7p 
describes a retarded propagator. It will be shown below 
that also r(g, qo;p,Po) can only be non-analytic in quad- 
rant II above the branch cut and IV below the branch 
cut in po and qq. We conclude that not only is it possi- 
ble to rotate the external energies, but we also move the 
contour of integration away from any singularities. This 
is illustrated in Fig. 1151 The only remaining singularities 
close to the integration contour are from the poles of the 
bosonic propagator as g ^ and A; ^ 0. This singularity 
is integrable and we will explicitly treat this below by a 
change of variables in Eq. (|C12[) . 

In the integral equation p3|) it is possible to integrate 
over the loop energy since Ek{q,pi,P2) is analytic in the 
lower half planes of both fermion energies. To see this, 
note that for any diagram contributing to S the two out- 
going fermionic propagators do not originate from the 
same bosonic propagator since then S would be two bo- 
son reducible. Instead each of the fermions originate 
from a boson which in turn decays into one other fermion 
which through a series of propagations forward in time 
contribute to the creation of the outgoing boson. Thus S 



contains only retarded propagators of —pi and —p2 and 
must be analytic in the lower half planes of the corre- 
sponding energies. 

It is not as simple to integrate out the loop energy in 
Eq. since Sfc((j, g',p— q') is not analytic in either 

(?g half planes. However, if one considers the integral 
equation satisfied by Sfc(q, q' ,p— q'), Eq. ([43]) with pi — 
q' and p2 = p — q' , the first product of propagators can 
be split into a term analytic in the upper half plane 
and a term analytic in the lower q'g half plane by noting 
that 

Go{K/2 + q - q')Go{K/2 -q-p + q') 
Go{K/2 + q-q')+Go{K/2-q-p + q') 
§-p,-{q--q')2/2-{q'+p-qr/2 + l0' 

The second product of fermionic propagators is treated 
similarly. The remaining part of the integral equation for 
Sfe(g, q' ,p — q') consists of a term analytic in the upper gg 
half plane and a term analytic in the lower q^ half plane. 
Thus effectively we have split S/c(g, q' ,p — q') into 

Ek{q,q',p-q') = E+ (q, q' ,p- q') + {q, q\p~ q') (C2) 

with {q,q' ,p — q') {E^ {q,q' ,p — q')) analytic in the up- 
per (lower) half plane. Note that this splitting means 
that in complete generality 

2fc(g,Pi,P2) S+(g,pi,p2) + 'E~{q,pi,p2). (C3) 

Sfc ('Zi'Z'iP ~ l') and E^{q,q\p — q') satisfy a set of 
coupled integral equations. The equation satisfied by 
(g,g',p- (?') is 



5fe {q,q',p~q') = 
An 
(2^) 



Go{K/2 + q~q') 



Go{K/2- 



f-po- (g - g')V2 - [q+p- 9*0 V2 + iO 
3 Go(i^/2 -Q- p)Do{K -Q- q')Ek{q, q' , Q) 



E 
2 



PO - {q + q'YI'i - {q- P + 9 ') V2 + M 



V2--E/4 



(C4) 



The right hand side of this equation is related to 
through Eq. (jCSp . The crucial point is that the 



right hand side of Eq. (jC4[) equals E^{q,p — q',q') 
which can be seen by writing the corresponding equa- 
tion for S^(g, q' ,p— q') and using the change of variables 
p - g' ^ g' along with Eq. (gS]). Thus S^(g,p - g', g') = 
"^k iq^q' 'P ^ q')- Consequently, using Eq. we con- 

clude that S^(g, g',p — g') and {q,q' ,p — q') has iden- 
tical contributions to F. We also note that 

Sfe(g,pi,P2) = 'B^{q,pi,p2) + Sj:(g,p2,Pi), (C5) 



which may be inserted in the right hand side of Eq. 
(|C4p . The result is an un-coupled integral equation for 
Sfc (g,g',p-g')- 



We now perform the gg-integration in the lower half 
plane in Eq. using Eq. (jCSp . This integration sets 
gg = g'^/2 — i?/4 > 0. At this value of the total outgoing 
fermion energy Eq. ()C4p reduces to 
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Sfe (g; q'., g"72 - E/A-p^ q',po - q'V2 + E / A) 



An 
1 

2^ 



GaiK/2 + q^q') 



Ga{K/2~q-q' 



Po + [q- 9*072 + {q+P~ q'Y/2 - E/2 - iO po + {q + q'? 12 + (-9 + p - g')V2 - S/2 - 
_^3^ S-(g; g", rV2 - £^/4; g/2 - i?/4) + (q; Q, Q/2 - i?/4; g", <f'V2 - £;/4) 



Po + QV2 + (P + Q) 72 - £^/2 - iO 1- A/-^ + 3Q74 + 3(f'74 + Q-g'/2 



(C6) 



The right hand side of this equation only contains 
evaluated at fcrmion energies related to the correspond- 
ing momentum by Qo — ~ E/A and q'^ = q'^/2 — E/A, 
which we will call "on shell". We will solve the integral 
equation by first letting 

Po^py2 + q'^ -p-q' -E/2, (C7) 

which means that the left hand side is also evaluated 
"on shell" . The resulting integral equation is an integral 
equation in \p — q'\, and the angle between these, 
using |fc|, IpI, \q\, and qo as input. Subsequently the "on 
shell" S~ can be inserted in the right hand side of Eq. 
(|C6p to evaluate S~ at any values of the outgoing energy. 

It is now possible to conclude that T(q, qQ]p,Pa) is ana- 
lytic in quadrant / and /// in po and go and that further- 
more any non-analyticity is far from the rotated contour 



of integration. The remaining pole of a fcrmionic prop- 
agator in Eq. after the go integration in the lower 
half plane is at po = 1 - P/A + g'2/2 + (p - q'f/2 - iQ 
where the real part is greater than 1 — fc^/4 > 1. In Eq. 
(IC6[) the poles in po are in quadrant // and also has real 
part less than — l-|-fc^/4 < —1. Finally, also in Eq. I|C6|) . 
the pole in go is in quadrant IV below the branch cut. 

We now define the "on shell" function 

2fe(g;pi,P2,cos6') = 

Sfc (g;pi,p-i 72 - E/A-P2J2 ' - E/A). (C8) 

Here, pi,2 = |pi,2| and 9 is the angle between these vec- 
tors. We then find the final set of equations to be solved 
to construct F: 



5fc(g;pi,P2,cos6') = - / ^ 

1 



q^+pl+pl-q-pi+q-p2-E-qo+pl/2 + {q-piY/2-E/2 

1 



g2 -Hp2 +p2 + g . Pi -q-p2-Eqo+ pl/2 + (g + pi)2/2 - E/2 
1 



d(cos0) / QMQ 



1 



1 - y/-E + ZQ'^/A + 3p74 + Qpi cos 0/2 
Sfc (g, pi , Q, cos 0) + Ek{q;Q,pi, cos 0) 



= .(C9) 



y/ [Q^ + Pi + P2 + QPi cos </) -I- P1P2 cos^ + Qp2 cos 6* cos — E)'^ — Q'^p^i^ — cos2 9)[1 — cos2 1 

We now insert Eq. I|C9|) in the right hand side of Eq. (|C6p using Eq. (jC8|) . and subsequently Eq. (|C6p in Eq. (Hi)) . 
Performing the integral over the azimuthal angle in the Q integral we find 



rfc(g,go;p,Po) = F[°^(g,go;p,po) + 



47^3 



d{cos 9) 



q''^dq' 



Po - g'2 - p2/2 pg' cos 9 + E/2 



X / d{cos(j)) 



Q^dQ 



1 - ^/-E + 3Q^/A + 3q'^/A + Qq' cos 0/2 
Sfc(g;g',Q,cos0) + Sfe(g; Q, g', cos0) 



7(Q2 -^p2/2 + Qp cos 6* cos -I- Po - E/2Y - Q^p^{l - cos2 9){1 - cos^ 0) 
where F^*^^ is the contribution from the first diagram on the right hand side in Fig. Illi which can be found to be 



(CIO) 



ri.*'^(g,go;p,Po) 



dq' 



1 



-g§ + {-E/2 + g'2 + p^/A + q^/A ~ q'pf 
167r2 Jo pg -E/2 + g'2 + p2/4 + g2/4 -gg + {-E/2 + g'2 -hp2/4 + q^/A + g'p)^ 



■In 
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X In ■ 



+ {-E/2 + q'^ + pV4 + qV4 - q'qy 



-pI + {-E/2 + 9'2 + p2/4 + gV4 + q'qf ■ 



(Cll) 



With the energies rotated to the imaginary axis all equa- 
tions used to find the t-matrix, Eq. (|42p at dq^ — > i dqo, 
and Eqs. (|C9IICll|) are real. This is not immediately 
obvious but may be easily checked numerically. 

To treat the poles in Eq. (|42)) of the bosonic propaga- 
tors at go = i(9^/4 — — iO) as k,q ^ we perform 
the following change of variables on the external energies 
and momenta to "polar" coordinates 



gV4 = Rl cos 01, 
pV4 = i?icos6l2, 



qo = Rf sin 6*1, 
Po^Rlsin02, (C12) 



where Ri £ [0,oo[ and 9i g [0, 7r/2]. Here we restrict the 
integration over q^ to the upper imaginary axis, using the 
symmetry of the integrand in Eq. (|42p as go ~^ ~Qo- This 
symmetry follows from the symmetry of the problem. For 
convergence reasons it is advantageous to integrate over 



finite intervals, and to this end we use the change of vari- 
ables 



Ri 



Zi + 1 



(C13) 



We employ the same change of variables on the internal 
loop momenta Q and q' in Eqs. (|C9|) and (|C10p . 



^2 ,2 

Q = 1, q = 1. 

^ ZQ + 1 ' ^ z' + l 



(C14) 



To solve the integral equations (|42|) and (|C9p we use the 
Nystrom method, writing the integral equations as ma- 
trix equations and inverting these to find the solution. 
Evaluation of the integrals in these integral equations is 
performed using Gauss-Legendre quadrature [27| . 
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